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Abstract 



We carried out a critical appraisal of the two theoretical models, Kurucz' ATLAS9 and 
PHOENIX/NextGen, for stellar atmosphere synthesis. Our tests relied on the theoretical fit 
^^ I of spectral energy distributions (SED) for a sample of 334 target stars along the whole spectral- 

f^ ■ type sequence, from the classical optical catalogs of Gunn & Stryker (1983) and Jacoby et al. 

(1984). The best-fitting physical parameters (Tcff, logg) of stars allowed an independent calibra- 

O ■ tion of the temperature and bolometric scale vs. empirical classification parameters (i.e. spectral 

type and MK luminosity class); in addition, the comparison of the synthetic templates from the 
ATLAS and NextGen grids allowed us to probe the capability of the models to match spectropho- 

^ I tometric properties of real stars and assess the impact of the different input physics. We can 

d • sketch the following main conclusions of our analysis: 

i) fitting accuracy of both theoretical libraries drastically degrades at low Toff, where both 
ATLAS and NextGen models still fail to properly account for the contribution of molecular 
features in the observed SED of K-M stars. 

t^ . ii) Comparing with empirical calibrations, both ATLAS and NextGen fits tend, in average, 

to predict slightly warmer (by 4-8%) Tcff for both giant and dwarf stars of fixed spectral type, 
but ATLAS provides in general a sensibly better fit (a factor of two lower a of fiux residuals) 
than NextGen. 

Hi) There is a striking tendency of NextGen to label target stars with an effective temperature 
and surface gravity in excess with respect to ATLAS. The effect is especially evident for MK I-III 
objects, where a fraction of stars of about one in four is clearly misclassified by NextGen in 
log (7. This is a consequence of some "degeneracy" in the solution space, partly induced by the 
different input physics and geometry constraints in the computation of the integrated emerging 
fiux (ATLAS model atmospheres assume standard plane-parallel layers, while NextGen adopts, 
for low-gravity stars, a spherical-shell geometry). A different T{t) vertical structure of stellar 
atmosphere seems also required for NextGen synthetic SEDs in order to better account for limb- 
darkening effects in cool stars, as supported by the recent observations of the EROS BLG2000-5 
microlensing event. 

Subject headings: stars: atmospheres - fundamental parameters 



1. Introduction 

Theoretical computation of model atmospheres 
has been a leading issue of stellar astrophysics 
in the last decades. In this framework, Kurucz' 
(1970; 1979) pioneering work certainly stands as 
a main reference, together with a few other ma- 
jor contributions like those of Gustafsson et al. 
(1975) and Tsuji (1976), on the synthesis of red 
giant stars. 

In its more recent versions, Kurucz' (1992a; 
1995) ATLAS code included over 58 million spec- 
tral lines, providing an accurate description of 
blanketing effects, that modulate ultraviolet emis- 
sion of stars (Holweger 1970; Gustafsson et al. 
1975), and also including nearly all the most im- 
portant di-atomic molecules, that shape spectral 
energy distribution (SED) at longer wavelength. 
The lack of tri-atomic molecules (in primis H2O), 
and an incomplete treatment of TiO opacity, how- 
ever, still prevents a satisfactory match of stars 
cooler than 3500 K (Kurucz 1992b; Castelh, Graf- 
ton, & Kurucz 1997). This limit of ATLAS the- 
oretical atmospheres unfortunately affects a num- 
ber of physical applications dealing, for instance, 
with the study of cool pulsating variables, or the 
match to integrated SED of galaxies, through stel- 
lar population synthesis. 

More recently, Hauschildt et al. (1999a, b) have 
presented their PHOENIX/NextGen grid of model 
atmospheres for dwarf and giant stars. AUard, 
Hauschildt, & Schweitzer (2000) extended the 
original bulk of models to the pre-main-sequence 
(pre-MS) evolution at the low-mass regime. With 
its 500 million atomic and molecular lines and 
a spherical geometry treatment of stellar struc- 
ture, the NextGen library is arguably the most 
advanced one currently available in the literature. 
As the low-temperature physical regime is suitably 
sampled, with models as cool as Toff — 2000 K, 
these may possibly fill the gap assuring a homoge- 
nous coverage of the stellar fundamental parame- 
ters across the whole H-R diagram. 

Given the relevance of the Kurucz and Hauschildt 
et al. contributions, it could be of special interest, 
at this stage, to carry out a combined analysis of 
the ATLAS vs. NextGen codes, in order to check 
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their mutual capabilities to match spectropho- 
tometric properties of real stars and assess self- 
consistency in their input physics. Our analysis 
follows the Hauschildt et al. (1999a) preliminary 
discussion, and will be carried out in two steps. 
After a brief description of the main features of 
each theoretical dataset (Sec. 2), we will first try 
a fit to observations of template stars and com- 
pare model outputs (Sec. 3). This will rely on 
an original optimization procedure, that also pro- 
vides an estimate of the fit uncertainty across the 
(Tefi, logg, [M/H]) phase space. The best fits will 
allow to establish an effective temperature scale 
and a calibration of the bolometric correction scale 
(Sec. 4). In Sec. 5 we will then analyse ATLAS 
vs. NextGen theoretical models for fixed fiducial 
spectral types. A full summary of the relevant 
conclusions of our tests is finally given in Sec. 6. 

2. Input model atmospheres and grid 
properties 

ATLAS model atmospheres assume steady- 
state plane-parallel layers under the hypothe- 
sis of local thermodynamic equilibrium (LTE). 
Line blanketing is computed statistically by 
means of opacity distribution functions (ODF), 
that average the contribution of the different 
atomic/molecular species through the correspond- 
ing oscillator forces (Strom & Kurucz 1966; Ku- 
rucz 1970, 1979). For our work, we used the AT- 
LAS 9 version (Kurucz 1995), whose treatment of 
convection is based on the mixing length theory 
(Bohm-Vitense 1958) and accounts for the so- 
called "approximate overshooting" , according to 
Castelli, Gratton, & Kurucz (1997). The mixing 
length parameter is set to £/Hp = 1.25 and the 
microturbulence velocity to 2 km s~^ throughout. 
The whole theoretical library is made available at 
the Kurucz web site.^ 

The model library spans a temperature range 
between 3500 < Tcff < 50 000 K, sampled at a 
variable step from 250 to 2500 K with increas- 
ing temperature; surface gravity and metallicity 
cover the interval 0.0 < logg < 5.0 dex at steps of 
Alog.g ^ 0.5 dex, and -5.0 < [M/H] < -hl.O dex, 
respectively. The corresponding SEDs, which also 
account for line opacity through the ODFs, span 
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from the far-ultraviolet (A — 90 A) to the far- 
infrared (160 /xm), sampled by 1221 wavelength 
points, with AA = 10 A in the UV and 20 A in 
the visual range. 

NextGen models have originally been computed 
by Hauschildt and collaborators with the multi- 
purpose code PHOENIX (e.g., Hauschildt ct al. 
1996). They assume LTE and plane-parallel ge- 
ometry for dwarf stars, while a spherical sym- 
metry is adopted in low-gravity model atmo- 
spheres (logg < 3.5) for giant and pre-MS stars 
(Hauschildt et al. 1999b; AUard, Hauschildt, & 
Schweitzer 2000). As a striking difference with 
respect to the Kurucz models, direct opacity sam- 
pling is performed including over 500 million lines 
of atomic and molecular species along the spec- 
trum. 

The phase-space domain of the NextGen grid 
spans the 2000 < T^s < 10 000 K range at steps 
ATcff = 100-200 K, with gravity in the interval 
0.0 < log (7 < 5.5 (Alogg = 0.5 dex), and metal- 
licity in the range -4.0 < [M/H] < +0.3. The 
SEDs cover the wavelength range from 100 A to 
970 /Ltm sampled at coarser wavelength steps, that 
can change from model to model (depending on 
the intervening absorption features in the spec- 
trum). However, a typical AA = 2 A step in the 
optical region can be picked out. All the data are 
available via anonymous ftp."^ 

In addition to the standard fundamental pa- 
rameters (i.e. Teff, \ogg, [M/H]), spherical models 
in NextGen require one supplementary "dimen- 
sion" in the phase space. As ^f oc M/B? and 
L oc R^T'^, then the emerging luminosity becomes 
L oc (M T'^)/g. Contrary to the plane-parallel 
case, therefore, mean surface brightness depends 
on the absolute size of stars through the stellar 
mass (M) . Giant-star models in the NextGen grid 
are computed for 5 Mq, while pre-MS stars as- 
sume 0.1 Mq. However, mass scaling is found 
to induce second-order effects on the stellar SED 
(Hauschildt et al. 1999b), although total luminos- 
ity of the models scales of course as L oc M for 
fixed temperature and gravity. 



^ ftp : //calvin . physast . uga . edu/pub/ 
http : //dllbert . physast . uga . edu/^yet i 
Note that the libraries of dwarf and giant stars available at 
these sites have lower T(,ff limits than the published ones. 



3. Matching SED of template stars 

A first "sanity" check in our analysis concerns 
the match with template stars along the whole 
O ^ M spectral-type sequence. Comparison with 
observed SEDs is one of the most natural applica- 
tion of model atmospheres. High-resolution spec- 
tra help in fact investigate chemical composition 
of stars, while the shape of (pseudo-)continuum at 
lower resolution gives clues on stellar gravity and 
effective temperature. 



3.1. The empirical spectral libraries 

For our test we considered the complete li- 
braries of stellar spectra by Gunn & Stryker (1983, 
hereafter GS83) and Jacoby et al. (1984, JHC84). 
Both sets of spectra span the whole range of stel- 
lar parameters, from giants (MK class I-III) to 
dwarfs (MK class IV-V), and have been widely 
used in the literature, particularly for population 
synthesis studies (e.g.. Pickles 1985; Guiderdoni 
& Rocca-Volmerange 1987; Fanelli, O'Connell, & 
Thuan 1987; Bruzual & Chariot 1993). The com- 
plete sample of target stars amounts to 336 objects 
(175 stars from GS83 and 161 from JHC84, with 
no stars in common to the two libraries). 

The GS83 data cover a wide wavelength range, 
from 3130 to 10800 A, observed at low resolution 
(FWHM=20 A in the blue and 40 A in the red) 
and sampled at steps of 10-20 A. Due to a poorer 
S/N quality in the ultraviolet, particularly for cool 
stars (Gunn & Stryker 1983), only the wavelength 
interval for A > 3500 A is suitable for our analy- 
sis. Thanks to a better sampling (AA = 1.4 A), 
the FWHM ~ 4.5 A resolution of JHC84 spec- 
tra is better exploited, giving a more detailed pic- 
ture of the main absorption features of template 
stars in the 3510 < A < 7427 A spectral range. 
For both the GS83 and JHC84 datasets we had 
to reject several regions (namely, around 6840- 
7000 A, 7140-7350 A, 7560-7720 A, 8110-8360 A, 
and 8900-9800 A) affected by telluric bands of O2 
and H2O. 

In our work we adopted the original MK spec- 
tral classification by GS83 and JHC84. For 24 out 
of 26 unclassified objects we relied on the SIMBAD 
database. All the spectra have been corrected for 
Galaxy reddening and atmospheric extinction as 
reported in the original data sources. 



3.2. The theoretical spectral libraries 

The subsample of [M/H] = model atmo- 
spheres has been used to match the observations. 
This choice is consistent with the mean metallic- 
ity of the GS83 and JHC84 stars. A systematic 
search from high-resolution abundance studies in 
the literature actually provided a mean value of 
[Fe/H] = -0.10 ± 0.24 for 67 stars in the GS83 
sample, and [Fe/H] = -0.06 ±0.19 for 25 JHC84 
stars.** 

A total of 409 theoretical SEDs have been 
collected from the ATLAS grid, while from the 
NextGen library we extracted a set of 314 theo- 
retical flux distributions composed by 164 plane- 
parallel models for dwarfs with 3000 < Toff < 
10 000 K and \ogg > 3.5 plus 150 spherical mod- 
els for giants with 3000 < T^s < 6800 K and 
logs < 3.5. The (Tcff, log5)-space coverage of 
the two grids is shown in Fig. 1. 

In order to consistently compare empirical and 
theoretical libraries, we degraded both GS83 and 
JHC84 spectra with a Gaussian kernel of FWHM 
= 25 A, rcbinning the output at constant steps 
of AA = 5 A. The same procedure has been ap- 
plied to the ATLAS and NextGen SEDs, sampled 
at the same set of wavelength points. The effect of 
this low-resolution approach on the results of our 
analysis (especially on the calibration of the tem- 
perature scale), will be discussed in some detail in 
Sec. 4 and 5. 



*The bulk of metallicity estimates for our star sample 
comes from the catalogs of Cayrel de Strobel et al. (1997, 
2001); other references are Hartkopf & Yoss (1982); Luck 
(1982); Kjaergaard (1984); Faber et al. (1985); Lambert 
et al. (1986); Norris (1986); Burkhart & Coupry (1988); 
Laird,Carney,& Latham (1988); Knude (1989); Luck & 
Bond (1989); Eggen (1991, 1998); Geisler, Claria, & 
Minniti (1991); Taylor (1991, 1999); Xu (1991); Friel & 
Janes (1993); Thogersen, Friel, & Fallon (1993); Worthey 
et al. (1994); Bartkevicius & Lazauskaite (1996, 1997); 
Claria, Piatti, & Osborn (1996); Zakhozhaj & Shaparenko 
(1996); Fry & Carney (1997); Flynn & Morell (1997); 
Schiavon, Barbuy, & Singh (1997); Takeda & Takada- 
Hidai (1998); Adclman (1999); Bartasiute, Ezhkova, & 
Lazauskaite (1999); Cenarro et al. (2001); Gray, Graham, 
& Hoyt (2001); Haywood (2001); Andricvsky et al. (2002); 
Venn et al. (2002), based on different spectroscopic or pho- 
tometric methods. Note that the only two metal-poor stars 
(HD 94028 and SAO 102986) in the JHC84 catalog have 
been excluded in our analysis so that we are eventually left 
with 334 objects in the GS83 plus JHC84 total sample. 
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Fig. 1.— The grid of 409 ATLAS and 314 
NextGen theoretical models considered in our 
analysis. A solar metallicity has been assumed 
throughout. 

3.3. Fitting procedure 

A "best fit" for the (Toff, logg) fundamental 
parameters (assuming [M/H] — 0) was searched 
for each star in the GS83 and JHC84 samples by 
matching the observed SED with both ATLAS and 
NextGen libraries. As described in full detail in 
Bcrtone (2001) and Bertone & Buzzoni (2001), our 
method basically relies on a minimization of the 
statistical variance in the relative flux domain, as a 
measure of the similarity between target spectrum 
and theoretical SEDs across the reference grid. 

Operationally, the spectrum of the i-th target 
star is compared with the j-th synthetic SED along 
the common wavelength range, deriving a residual 
function 

X(,,,)(A)=ln/,(A)-ln/,(A)+k (1) 

in the flux logarithm domain, as shown in Fig. 2. 
The offset constant, k, in eq. (1) is such as 



^X(,,,)(A)^0, 



(2) 



so that 

k=<[ln/,(A)-ln/,(A)]>, (3) 

while the standard deviation 

s(X)(,,,) = VVar[X{\)] (4) 
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Fig. 2.— The observed SED of star no. 45 from GS83 (namely, HD 154760, of spectral type G2V) is 
compared with several ATLAS models. Left plots are for fixed gravity and AToff — ±250 K around the 
reference value of the central panel. Right plots are instead for fixed Toff and Alogg = ±0.5 dex. At 
the bottom of each panel we display the residual function X(^ij-^{X) according to eq. (1) with its standard 
deviation, s, from eq. (4) as labelled. The central panel is the ATLAS best-fit solution for this star. 
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Fig. 3. — Standard deviation of flux residuals for NextGen (left panel) and ATLAS (right panel) fitting 
models for star BD +38 2457. Each curve connects equi-gravity points in the grid vs. Tpg (dotted line for 
log (7=0, dash-dotted for logg^c ~ 5.5 or log^ATLAs = 5.0, dashed for intermediate values). Solid line is 
the best-gravity solution with its minimum marked by the big open square. The horizontal dotted line 
shows the upper limit for Smin at a 2-g confidence level, as resulting from an F statistical test. The fiducial 
fundamental parameters for this case are T^fi— 4850]']™q K, logg=5.5 ± 0.5 dex with the NextGen models. 
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K and log5=4.0±0.5 dex with the ATLAS models. 
0.06) compared to NextGen (smin = 0-12). 



Note the better accuracy of the 



provides a measure of the spectral likelihood be- 
tween observations and theoretical models.^ The 
underlying hypothesis of our approach is of course 
that a non-degenerate trend exists for the s(A') 
function in the (Teff, log.g) phase space, so that a 
univocal best solution can be found for a given in- 
put SED. For each target star we then mapped the 
s(X) distribution by matching the whole grid of 
synthetic model atmospheres, and searched for an 
absolute minimum, Smin, after performing a cu- 
bic spline interpolation of the s{X) points at each 
gravity. This allowed us to locate the best-fitting 
values of Teff and log g with a nominal resolution of 
ATeff ~ 10 K and Alogg ^ 0.5 dex, respectively. 
Statistical uncertainty of fiducial distinctive pa- 
rameters was estimated by means of a one-tail F 
test on the value of Smin, at a 95% confidence 
level. ^ The Smin confidence interval was trans- 
lated into an equivalent (ATeff, Alogg) error box 
relying on a first-order estimate of ds{X)/dTes 
and ds{X)/ dlogg evaluated around the Smin re- 
gion in the phase space. An example of the fitting 



®The freedom degrees for s(X) are settled by the number 
A'^ of wavelength points after spectrum broadening, as de- 
scribed in previous section. 

^The freedom degrees of the F distribution, in this case are 
simply N — 1, cf. previous footnote. 



procedure for a star in the GS83 sample is dis- 
played in Fig. 3. 

The robustness of our minimization procedure 
was probed by means of a bootstrap test. We 
added a 10% noise to the full set of Kurucz syn- 
thetic SEDs and tried our best fit to recover the 
original (i.e. unperturbed) reference parameters. 
In all cases, the correct Toff was identified, typi- 
cally with a 1-2% uncertainty, a value that raised 
to 5-8% just for the few poorest cases. As far as 
surface gravity is concerned, the nominal values 
of the reference models were picked up within a 
±0.5 dex in 98% of the cases (i.e. with only 8 
outliers out of 409 fitted SEDs). 

4. Temperature scale calibration 

Out of the total of 334 stars in the GS83 and 
JHC84 libraries, a consistent fitting solution for 
the (Toff', log g) fundamental parameters was found 
for 272 and 230 stars, respectively, using ATLAS 
and NextGen reference grids. Most of the remain- 
ing unfitted objects are M and 0-B stars, that 
is at the two extreme edges of the temperature 
scale, where a fair value for Smin cannot be con- 
fidently located within the theoretical model grid. 
The NextGen code, however, has proven to be 
marginally more efficient in the fit of M stars (of a 
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Fig. 4. — The Smin distribution vs. spectral type 
for the stars of the GS83 (upper panel) and JHC84 
atlas (lower panel) . Open markers indicate the re- 
sults for ATLAS models, filled symbols those from 
the NextGen grid. Note, for the latter, that a limit 
at Teft< 10 000 K does not allow any fit to 0-B 
stars. The value of Smin is a measure of the mean 
percent accuracy of the best fit to the observed 
spectrum. 



total of 38 stars in this class, 10 were successfully 
matched by ATLAS and 24 by NextGen). 

The accuracy of the ATLAS and NextGen 
model libraries in the fit of GS83 and JHC84 stars 
can be analysed by means of Fig. 4. In the two 
panels of the figure we report the distribution of 
the residual standard deviation of the best fits for 
stars in the two observed samples. As a common 
feature in the two plots, note that NextGen pro- 
vides in general a poorer fit compared to the AT- 
LAS code. This is particularly evident for G and 
K stars (Tes ~ 5500 -^ 4000 K), where ATLAS is 
a factor of two better than NextGen in terms of 
best-fitting variance. The figure also shows that 
the accuracy in the definition of the temperature 
scale directly depends on the wavelength baseline 
of the spectra. Compared with the GS83 stars, 
in fact the JHC84 fits are slightly poorer, given 
a narrower spectral range for the JHC84 library 
(i.e. AA ~ 4000 A vs. 7500 A for GS83). 

The incomplete treatment of molecular opacity 
in the Kurucz code is well evidenced in both plots 



of Fig. 4, facing the sharp increase of standard 
deviation in the fit of late-K and M stars. The 
same problem seems to affect, at a similar level, 
also the NextGen fits confirming in any case a still 
unsolved and pervasive discrepancy of the theory 
to self-consistcntly reproduce cool stars. 

The temperature scale resulting from the AT- 
LAS and NextGen fits is displayed in Fig. 5. 
We compared with a number of empirical mean 
loci for dwarfs and giants including the John- 
son (1966) classical compilation and the Bohm- 
Vitense (1981) scales for hot stars. We also con- 
sidered the recent calibration of F0-K5 giant stars 
from Alonso, Arribas, & Martinez-Roger (1999) 
and the extension to late-M giants of Perrin et al. 
(1998), which includes the Ridgway et al. (1980) 
data. Furthermore, Di Benedetto (1998) provided 
accurate and systematic measures of the effective 
temperature for a wide sample of 537 stars of A 
to K spectral type, within an internal accuracy 
of ±1% in the individual Tcn estimates, using the 
surface brightness technique (Wesselink 1969), cal- 
ibrated by means of the angular diameters of 22 
stars. His mean locus for the dwarf and giant sub- 
samples is superposed to our data in Fig. 5. 

In general, we find a consistent trend between 
the Teff scale from the theoretical fits of the GS83 
and JHC84 stars and the empirical reference cal- 
ibrations. A more scattered distribution of our 
points derives, of course, from the fact that these 
are fits of individual stars, instead of a mean lo- 
cus, and because of the luminosity and metallicity 
spread of the sample. Only one clear outlier ap- 
pears among the GS83 stars in the top-right panel 
of Fig. 5; this is star HD 156252 (ahas 38 Oph), 
classified by GS83 as a type AlV star with (dered- 
dened) (B-V) = —0.14 mag and color excess E(B- 
V) = 0.16 mag. The exceedingly blue color calls 
for a warmer fitting temperature (Toff ~ 14 000 K) 
with ATLAS, while also the NextGen match sug- 
gests a temperature in excess of 10 000 K, as no 
minima of the standard deviation s [see eq. (4)] 
were present below that temperature. 

This star is reportedly among the most red- 
dened ones in the GS83 list but if one accounts, 
alternatively, for a much lower color excess as re- 
ported in the Hipparcos catalog [namely E(B-V) — 
0.015 mag], then the evident discrepancy between 
fitting temperature and spectral type could easily 
be recovered. 
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Fig. 5. — The effective temperature scale as derived from the GS83 (upper panels; dot markers) and JHC84 
stars (lower panels; diamond markers), after NextGen (left panels) and ATLAS (right panels) fits. Our 
results are compared with several empirical calibrations from the literature, as indicated in the top left 
panel. The outlier star HD 156252 is indicated in the top right panel (note that this object is not comprised 
in the NextGen plots due to an upper limit to effective temperature of the model grid at Toff < 10 000 K). 
See text for discussion. 



To better single out the differences between the 
NextGen and ATLAS fits, in Fig. 6 and 7 we con- 
sidered separately the main sequence (MK class 
V) and giant (MK III) star subsamples display- 
ing the ATeff /TeS between our fitting tempera- 
ture and the reference calibration of Bohm-Vitense 
(1981) and Johnson (1966), for dwarfs, and Alonso 



et al. (1999) and Perrin et al. (1998), for giants. 
Fig. 6 shows that both ATLAS and NextGen grids 
tend to fit F to M stars with a 4-8% warmer effec- 
tive temperature, the Teff excess being in general 
higher for the NextGen fits. 

The situation is somehow different for giants 
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Fig. 6. — Teff residuals vs. spectral type of ATLAS 
and NextGen best-fits for the subsample of MK V 
stars and the corresponding calibration of Johnson 
(1966), for A-M stars, and Bohm- Vitense (1981) 
for 0-B-types. Open dots identify GS83 objects, 
while diamonds mark the JHC84 stars. The out- 
lier HD 156252 is off the ATLAS plot, as indicated 
by the vertical arrow (see text for discussion). 



(see Fig. 7) with a drift in the point distribution 
with respect to the Alonso ct al. (1999) and Per- 
rin et al. (1998) Teff calibration (but, again, with 
the NextGen output marginally warmer than the 
ATLAS one). 

As for the Kurucz models, a glance to Fig. 2 
makes clear that part of the bias toward higher 
fitting temperatures might derive from the blan- 
keting effects in the ultraviolet region of the stellar 
SED, shortward of 4000 A. The residual scatter in 
this spectral region is in fact a major source to 
the global variance when matching observed SED 
and theoretical models, thus sensibly constraining 
the choice of the best-fit solution. More entan- 
gled is the situation for NextGen models, which 
adopt a different chemical mix representative of 
solar metallicity. While in fact ATLAS relics on 



Fig. 7. — As for Fig. 6, but for MK III giant stars 
vs. the corresponding calibrations of Alonso et al. 
(1999) for F-K stars, and Perrin et al. (1998) for 
M-type stars. 



the Anders & Grevesse (1989) solar abundances, 
NextGen assumes the revised values by Jaschek & 
Jaschek (1995); for Z — Zq, this makes NextGen 
Fe abundance slightly lower compared to ATLAS 
(namely, [FeNc/FeATLAs] - -0.17 dex). How- 
ever, this feature could hardly explain the ob- 
served trend in the Tcs calibration as a Fc-poorer 
model atmosphere should actually display a lower 
blanketing, thus allowing a cooler temperature to 
fit the observed SED of stars (see, in this sense, the 
quantitative discussion by Buzzoni ct al. 2001).^ 



^As a cross-check in this regard, we tried a fit of the GS83 
stars relying on the Kurucz Ubrary with [Fe/H] = —0.1 
instead of solar. From the operational point of view, this 
should roughly mimic the NextGen solar case. As expected, 
the GS83 fitting temperatures are in average 50—100 K 
cooler than the values obtained with the ATLAS library 
at [Fe/H] = 0. 
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Fig. 8. — Derived bolometric correction for GS83 
(open dots) and JHC84 (diamonds) stars accord- 
ing to ATLAS and NextGcn model atmospheres. 
Our results are compared with the empirical cali- 
brations of Flower (1977, dotted line for MS stars, 
dashed for giants), Johnson (1966, solid), Bessell 
(1991, dashed-dotted), and Code et al. (1976, cross 
symbols). 



4.1. Bolometric corrections 

Our best-fitting procedure with model atmo- 
spheres allows in principle a straightforward es- 
timate of bolometric luminosity for stars in the 
GS83 and JHC84 samples. When coupled with 
the individual V magnitudes, this could eventu- 
ally supply a measure of the bolometric correction 
(BC). 

In order to set the BC scale, for our calcula- 
tions we identified the theoretical template for the 
Sun as the ATLAS model SUNK94, with (Tes, 
log 5, [M/H]) == (5777 K, 4.44, 0.0) according to 
Castelli, Gratton, & Kurucz (1997). Its theoretical 
SED has been convolved with the V filter profile 
of Bessell (1990) and photometric zero points were 
tuned up such as to have BCq — BoIq — Vq = 
-0.07 (Bessell, Castelh, & Plez 1998). 

The BC for each star in the GS83 and JHC84 



samples is then computed as: 

BC = -2.5\og{aT;^fj) -V -k + 2.22. (5) 

Note that in eq. (5), the V magnitude derives from 
the convolution of the observed spectrum while 
the offset "k" (that properly scales the bolometric 
magnitude of the theoretical fitting SED) is from 
eq. (3). 

Our results are compared, in Fig. 8, with other 
standard calibrations for dwarf and giant stars 
vs. spectral type. We considered in particular 
the work of Flower (1977) and Johnson (1966) 
and its later revisions of Bessell (1991) (for late- 
K and M dwarf stars), and Code et al. (1976) 
(for hot 0-B stars). When necessary, bolometric 
scales were shifted consistently such as to assure 
BCq = -0.07. 

The final reference calibration for Toff and BC 
vs. spectral type for dwarf and giant stars ac- 
cording to ATLAS and NextGen is summarized 
in Table 1. Given the limited spectral coverage 
of the observed spectra, our bolometric extrapo- 
lation suffers, of course, from intrinsic uncertain- 
ties at the two extremes edges of the temperature 
scale, where a substantial fraction of stellar en- 
ergy is emitted outside the optical range. The BC 
calibration of Table 1 for B and M stars should 
therefore be taken with some caution since it crit- 
ically relies on the theoretical input physics. 

A more extensive discussion of this issue, from 
a fully theoretical point of view, has been carried 
out by Bessell, Castelli, & Plez (1998) based on 
the ATLAS and NMARCS model predictions, the 
latter ones as from the calibrations of Plez et al. 
(1992) for giants and Edvardsson et al. (1993) for 
dwarfs. A comparison of our results with those of 
Bessell, Castelh, & Plez (1998) is shown in Fig. 9 
comforting, however, on the general agreement of 
our output. 

5. Comparing template SED along the 
spectral-type sequence 

The whole set of synthetic templates for the 
GS83 and JHC84 stars prompts a straightforward 
comparison of the ATLAS vs. NextGen code per- 
formances taking into account in a self-consistent 
way the effect of the different input physics on the 
match of SED for real stars along the O -^ M 
spectral-type sequence. For our test we especially 
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Table 1 
Fiducial ATLAS and NextGen calibration for temperature scale and bolometric 

CORRECTION 
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NextGen 
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Giants 
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Giants 
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Tcff 


BC 
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BC 


Tcff 


BC 


Tcff 
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BO 
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0.05 
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0.02 
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GO 
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Fig. 9. — As for Fig. 8, but vs. temperature scale for several theoretical calibrations. Our results, from Table 
1 (solid lines), are compared with those of Bessell, Castelli, & Plez (1998), also based on the ATLAS models 
(open dots), Plez et al. (1992) for cool giant models (crosses), and Edvardsson et al. (1993) for M dwarfs 
(triangles), both relying on the NMARCS atmosphere code of Gustafsson et al. (2003). Left panel reports 
our NextGen calibration, while right panel is for the ATLAS theoretical locus. 



relied on the subset of 216 objects from the GS83 
and JHC84 catalogs for which a nominally best fit 
exists both for the ATLAS and NextGen grids over 
the temperature range 3500 < Teff < 10 000 K. 

A first interesting feature, when comparing the 
two sets of theoretical templates, like in Fig. 10, 
concerns the distribution of the fitting gravity. It 
is evident from the histogram that for an impor- 
tant fraction of target stars NextGen tends to fit 
with a higher gravity with respect to ATLAS. This 
actually led to a number of "catastrophic outliers" 
among the giant and supergiant stars in the GS83 
and JHC84 samples, as shown in Fig. 11. In the 
NextGen plot, in fact, (upper panel of the figure) 
23 out of the 100 MK I-III stars are unexpect- 
edly located in the high-gravity region of the dia- 
gram, pertinent to class V dwarfs, with a nominal 
"best-fit" gravity of \ogg — 5.5 dex. Conversely, 
only four such gravity outliers are present in the 
ATLAS diagram (lower panel) with a gravity of 
log 5 = 5 dex. 

A similar trend can also be recognized for the 
fitting stellar temperature T^g, as we have been 
discussing in Sec. 4. Again, Fig. 12 shows that 
NextGen Tcs estimates are in average 2% in ex- 



cess with respect to the ATLAS best-fit values, 
with a sensibly higher scatter for the JHC84 stars 
(see lower panel in the figure), partly depending on 
the shorter wavelength baseline compared to the 
GS83 set of spectra [cr(Arcff7T'cff) ^ 0.031 dex for 
the JHC84 sample vs. a value of 0.018 dex for the 
GS83 stars]. 

The tendency of NextGen to overestimate tem- 
perature and gravity can be illustrated by means 
of Fig. 13, where we map the distribution of the 
fit variance (s) across the theoretical grid for two 
MK III giants in the GS83 and JHC84 samples.*^ 
One sees from the plots that actually two phys- 
ically distinct solutions exist for these stars, one 
that correctly locates both K2III giants in the 
low-temperature low-gravity range [namely, (T^g, 
\ogg) ^ (4500 K, 3.0 dex) in our example] and 
adopts a spherical model (i.e. in the logg < 3.5 do- 
main), and the other (nominally better) one that 
assumes a plane-parallel geometry but places stars 
at a much higher logg ~ 5.5 dex and Ten ^ 
4800 K. This apparent "bimodality" in the solu- 



'This plot is basically a projected view of the 3-D fitting 
surface, like that shown in Fig. 3 
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Fig. 10. — The log^Nc vs. log^ATLAs distribution 
for the sample of 216 stars with common fitting 
solutions. The vertical axis shows the frequency 
number. Note the excess of high-gravity best- 
fit solutions for NextGen, compared to the cor- 
responding ATLAS distribution (sec text for full 
discussion). 



tion space is also well evident in the J11C84 panel 
of Fig. 12, where giant star distribution appears to 
split in two distinct sequences depending whether 
the spherical or the plane-parallel solution prevails 
as a best fit.*^ 

Note, by the way, that some correlation in the 
temperature and gravity excess, when fitting em- 
pirical SED with theoretical models, can natu- 
rally be expected on the basis of the arguments 
pointed out by Buzzoni et al. (2001). Their ex- 
periments showed in fact that a correspondingly 
higher gravity should likely be required to recover, 
at medium-high resolution, "too shallow" absorp- 
tion features as predicted by a too warmer model 
forced to match the low-resolution SED of a given 
star. According to Buzzoni ct al. (2001), such a 
tight dependence between A log Toff and Alogg 
can be written in the form 



^^^aOOofi^VdexK- 

AlogTeff \TeS J 



(6) 



o 3 
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M „ 

_o 3 
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ATLAS 




4.6 



4.4 



4.2 4.0 

log T,„ 



3.6 



^The same effect is not equally evident in the GS83 plot due 
to the dominant fraction of dwarf stars in this sample. 



Fig. 11. — The log Teff vs. log 5 distribution of 
ATLAS and NextGen fitting output for GS83 and 
JHC84 stars. The different symbols mark the lu- 
minosity class: MK V = •, IV = big •, III = o, 
II-I — Q. Overplotted are also the evolutionary 
tracks for stars of solar metallicity and M = 1,2, 
and 5 Mq from Girardi et al. (2000), and for 
10 Mq from Salasnich et al. (2000). The white 
region show the parameter space covered by the 
theoretical libraries (see also Fig. 1). 



A consistent trend in this sense is confirmed in 
Fig. 14, that reports the temperature and gravity 
differences between ATLAS and NextGen fiducial 
solutions for the 216 target stars in common. 

5.1. Sphericity effects on theoretical SED 

Such a different behaviour of ATLAS and 
NextGen model output directly calls, of course, 
for a distinct physical approach in the calculation 
of the inner structure of the stellar atmosphere. 
This is especially true for giant stars, where the 
plane-parallel model atmospheres of ATLAS are 
to be compared with the spherical-shell geome- 
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Fig. 12. — Temperature residuals for the 216 GS83 
(upper panel) and JHC84 stars (lower panel) with 
both ATLAS and NcxtGcn best-fit solution. Stars 
arc labelled according to their MK luminosity class 
(MK IV-V: open markers; MK I-III: solid mark- 
ers). Note, in the JHC84 plot, the peculiar distri- 
bution of giant stars, along two distinct point se- 
quences. Over the whole sample, NextGen tends 
to predict, in average, an effective temperature 
about 2% warmer than the ATLAS value (see text 
for a full discussion of these two important fea- 
tures). 



try of NextGen. The impact of geometry on the 
emerging flux of the theoretical models has been 
first assessed in a pioneering work by Scholz & 
Tsuji (1984) on the atmospheres of M and C stars, 
and more extensively explored in the recent years 
by Plez and collaborators (Plez 1990; Plez et al. 
1992). 

Basically two intervening effects modulate the 
integrated SED of spherical models with respect to 
their corresponding plane-parallel cases. First, as 
a general trend for fixed T^s and logy, spherical 
model atmospheres tend to display a lower elec- 
tronic pressure and a cooler temperature profile 
vs. stellar spatial coordinate (i.e. radius or opti- 
cal depth, cf. e.g. Scholz & Tsuji 1984). To some 
extent, this is the physical consequence of the de- 
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Fig. 13. — An illustrative example of the NextGen 
fitting procedure for two K2 giant stars from the 
GS83 and JHC84 samples. Plotted is the map 
of the standard deviation, s, of residual flux be- 
tween observed and theoretical SED across the 
model grid [according to eq. (4)] in the Tcff vs. logy 
phase space. One sees that two best-fit solutions 
can be identified in each plot, one correctly plac- 
ing the K2 III stars in the low-temperature low- 
gravity region (i.e. Tcs ~ 4500 K, logy ~ 3.0 dex) 
and the other nominally better one shifting stars 
to slightly warmer temperature and much higher 
gravity (Teff ~ 4800 K, logy ~ 5.5 dex). 



creasing gravity when moving outward of stellar 
photosphere; with a lower gravity, in fact, thermo- 
dynamical equilibrium in the external layers read- 
justs such as to allow a lower pressure of the elec- 
tronic plasma (because of an increased mean dis- 
tance between atoms and a higher dumping poten- 
tial for the bound-bound and bound-free e^ tran- 
sitions) and a cooler temperature, still sufficient 
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however to "sustain" the atmosphere structure. 

As a result, for fixed Tcs and logg, the SED of 
a "spherical" star is therefore expected to display 
sharper absorption lines and a "redder" contin- 
uum. Among others, this should also reflect in a 
less severe blanketing absorption (see in this sense 
the experiments of Hauschildt et al. 1999b) as a 
consequence of a reduced blend of metal absorp- 
tion lines at short wavelength. 

A second related effect that should be dealt 
with, when comparing plane-parallel and spher- 
ical model atmospheres, concerns limb darken- 
ing. Due to the geometry, in fact, the integrated 
flux that emerges from a "spherical" star receives 
a more important contribution from low-gravity 
cooler layers, and appears therefore in average 
"cooler" with respect to its corresponding plane- 
parallel model (Claret & Hauschildt 2003). 

This feature also emerges from the recent re- 
sults of Fields et al. (2003) on the microlensing 
surface scanning of the K3 giant star related to 
the EROS BLG2000-5 event. Surface brightness 
measurements for this star are in fact inconsistent 
with the NextGen best-fit predictions at higher 
than 10(7, and indicate that the derived T{t) ver- 
tical structure of the theoretical atmosphere sen- 
sibly overestimates limb-darkening effects. 

6. Summary and conclusions. 

In this work we carried out a combined compar- 
ison of the two theoretical codes ATLAS (Kurucz 
1992b) and NextGen (Hauschildt et al. 1999a,b) 
for stellar atmosphere synthesis. Our tests relied 
on the fit of a set 334 target stars of nearly solar 
metallicity, spanning the whole sequence of spec- 
tral types and luminosity class, observed in the 
optical range by Gunn & Stryker (1983) and Ja- 
coby et al. (1984). 

For about 80% of this sample we obtained an 
estimate of the physical parameters {Tcs, logg) 
of stars and their related statistical uncertainty 
by means of an original fitting procedure that 
matched the observed SED with the ATLAS and 
NextGen model grids. 

This provided a twofold application of our re- 
sults; from one hand, we achieved a self-consistent 
and independent calibration of the temperature 
and bolometric scale for giant and dwarf stars 
vs. empirical classification parameters (i.e. spec- 




-0.02 



-0.01 0.00 0.01 0.02 0.03 0.04 
log T,„,NG - log T,„.4Tus 



Fig. 14. — Temperature and gravity difference be- 
tween ATLAS and NextGen best-fit solutions for 
216 stars in the GS83 and JHC84 samples. Marker 
size is proportional to the MK luminosity class (i.e. 
big markers — MK I-III giants, small markers = 
MK IV-V dwarfs). It is evident a correlation be- 
tween ATeff and A log 17, especially for giant stars 
in the JHC84 plot. 



tral type and MK luminosity class). On the 
other hand, the comparison of the synthetic tem- 
plates from the ATLAS and NextGen model grids 
allowed us to directly assess the relative per- 
formances of each theoretical code to reproduce 
SED of real stars according to the different input 
physics adopted. 

The comparison of our results with several em- 
pirical calibration scales in the literature (see Sec. 
4) led to the following main conclusions: 

i) the good fitting accuracy (ctaux ^ 2 — 5%) 
of both theoretical models in reproducing SED 
of early-type stars (spectral type F and earlier) 
drastically degrades at lower Toff, especially for K 
stars, where both ATLAS and NextGen codes still 
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fail to properly account for the increasing contri- 
bution of molecular features in the spectra of stars. 
In general, ATLAS is found to provide a systemat- 
ically better fit (a factor of two lower residual ctaux) 
than NextGen along the whole B ^ K spectral- 
type sequence, although the NextGen grid, due 
to its lower Toff limit, more efficiently matches M 
stars. 

ii) Comparing with empirical calibrations, both 
ATLAS and NextGen tend, in average, to pre- 
dict warmer (by 4-8%) T^g for both giant and 
dwarf stars of fixed spectral type. As for the AT- 
LAS models, this effect has probably much to do 
with the imperfect treatment of metal blanketing 
at short wavelength (as extensively discussed, for 
instance, by Castelli, Gratton, & Kurucz 1997), 
while the case of NextGen seems more entangled. 

This issue has been further explored in Sec. 5 
by comparing the ATLAS vs. NextGen template 
sequences for 216 stars in the GS83 and JHC84 
catalogs with nominal best fit in the 3500 < T^g < 
10 000 K temperature range. As a general feature, 
NextGen best-fit solutions are found to predict a 
temperature and gravity excess with respect to the 
corresponding ATLAS solutions for given target 
stars. The effect is especially evident for MK I-III 
objects, where the NextGen fails to correctly set- 
tle log g and sensibly overestimates surface gravity 
in about 25% of the cases vs. 4% of ATLAS. This 
misclassification partly derives from a lower capa- 
bility of NextGen spherical models to reproduce 
the SED of giant stars, compared to the fit with 
the plane-parallel geometry. In most cases the lat- 
ter proved in fact to be formally more accurate, 
leading however to a less physical combination of 
the fundamental parameters of stars. 

An in-depth analysis of the fit accuracy for SED 
of target stars shows that, to some extent, the 
NextGen Tes and \ogg excess is correlated, as a 
consequence of a sort of "degeneracy" in the so- 
lution space (Buzzoni et al. 2001). The effect is 
likely magnified in our framework when consider- 
ing that, for low-gravity stars, ATLAS model at- 
mospheres assume standard plane-parallel layers 
while NextGen adopts a spherical-shell geometry. 
Because of a more important contribution of exter- 
nal atmosphere layers to the integrated emerging 
flux in the case of the NextGen output, for fixed 
Teff and logc; this implicitly calls for a "redder" 
theoretical SED and a reduced blanketing absorp- 



tion of metal blends at short wavelength, as a con- 
sequence of "sharper" spectral features. 

The possible overestimate of the limb-darkening 
effects, as a consequence of the adopted T{t) ver- 
tical structure of NextGen model atmospheres, 
seems also a critical issue in this regard, as in- 
dicated by the recent observations of the EROS 
BLG2000-5 microlensing event. 
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